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Abstract We report on initial results of our Spitzer Cycle 2 
program to observe the young massive star R Mon and its 
associated HH 39 Herbig-Haro object in the mid-infrared. 
Our program used all instruments on-board Spitzer to ob- 
tain deep images with IRAC of the HH 39 complex and of 
R Mon and its surroundings, a deep image of HH 39 at 24 
and 70 jim with MIPS, and mid-infrared spectra with the 
SH, LH, and LL modules of IRS. The aim of this program is 
to study the physical links in a young massive star between 
accretion disk, outflows and jets, and shocks in the associ- 
ated HH object. Our preliminary analysis reveals that several 
knots of HH 39 are clearly detected in most IRAC bands. 
In IRAC4 (8 jim), diffuse emission, probably from PAHs, 
appears as foreground emission covering the HH 39 emis- 
sion. The HH 39 knots are detected at 24 microns, despite 
the fact that dust continuum emission covers the knots and 
shows the same structure as observed with IRAC4. The IRS 
spectra of HH 39 show weak evidence of [Ne II] 12.8 jim 
and 0-0 S(l) H2 17.0 jim lines. A more detailed analysis is, 
however, required due to the faintness of the Herbig-Haro 
knots. Finally, we obtained the SH and MIPS SED spectra 
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of R Mon. A PAH emission feature at 11.3 jim is detected 
on top of the strong continuum; although no strong emission 
or absorption lines are observed, we will seek to detect faint 
lines. The combined IRAC, IRS, and MIPS data of the R 
Mon/HH 39 system will help us to understand circumstellar 
disk processing, and the connection between jets, outflows, 
and HH objects. 

Keywords Herbig Be star • Herbig-Haro Object • Accre- 
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1 Introduction 
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Herbig Ae/Be stars (hereafter HAEBEs; Herbig, 1960) form 
a class of massive, young stars with high luminosities (10 — 
1000 Lq) and with a strong infrared (IR) excess due to cir- 
cumstellar dust (Waters and Waelkens, 1998). The spectral 
energy distribution (SED) of HAEBEs can be explained by 
circumstellar envelopes with polycyclic aromatic hydrocar- 
bons (PAHs; 3.3, 6.2, 7.7, and 1 1.3 ^m; Brooke et al, 1993; 
Meeus et al., 2001), and by dust emission from amorphous 
or crystalline silicate bands (8 — 12 /im) or molecular and 
atomic transitions. Accretion disks and outflows in the less 
massive classical T Tau stars (CTTS) are intimately con- 
nected, thus the presence of outflows in some HAEBEs sug- 
gests that disks could be present as well (e.g., Corcoran and 
Ray, 1998a). Corcoran and Ray (1998b) also found that the 
wind mass-loss rate correlates with the IR excess over 5 
orders of magnitude in luminosity and from 0.5 to 10 M 
when using both CTTS and HAEBEs. Outflows in HAEBEs 
are ~ 2 — 3 times faster (v ~ 600 — 900 km s _1 ) than in 
CTTS, but generally show similar collimation (3° — 10°), al- 
though the fraction of poorly collimated outflows in HAEBEs 
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(50° - 120°) is larger (e.g., Mundt and Ray, 1994). Complex 
shocks occur at the interface between the jet and the molecu- 
lar material (Draine, 1980), heating up the gas which in turn 
cools down radiatively. This gas is detected as Herbig-Haro 
objects (HH) in excited lines in the optical, in the near-IR 
(e.g., H2 ro- vibrational lines around 2 /im), and in the mid- 
IR ([O I] 63 Aim; Nisini et al, 1997; Liseau et al., 1997; 
Molinari et al, 2000). 

2 The R Mon and HH 39 System 

The Herbig Be star R Mon (d — 800 pc) is associated with 
NGC 2261, a reflection nebula that gradually faints with in- 
creasing wavelength (Close et al., 1997). A bipolar outflow 
and a high- velocity jet (v ~ 100 km s _1 ) pointing toward 
the nearby HH 39 knots (r = 7.5', PA = 350°) were detected 
(Canto et al., 1981; Brugel et al., 1984; Movsessian et al., 
2002). Evidence for a circumstellar disk is substantiated as 
well (e.g., Beckwith et al., 1986; Fuente et al, 2003; Fuente 
et al., 2006). A faint CTTS companion at a separation of 
0.7" and with K' flux of 34 mJy is mentioned by Close et 
al. (1997) and Weigelt et al. (2002). The HH 39 cluster of 
knots covers an area about 25" x 45" (Fig.[TJ Jones and Her- 
big, 1982; Walsh and Malin, 1985), and the northernmost 
knot, knot A, is believed to be the working surface of the jet 
(Schwartz and Schultz, 1992). 

Near-IR spectra of R Mon exist mostly in the 0.9 — 2.5 /im 
range with low (~ 400) spectral power (Kelly et al., 1994; 
Reipurth and Aspin, 1997; Porter et al., 1998). ISO neither 
observed R Mon nor HH 39. The IRAS LRS spectrum (9 - 
21 jUm) of R Mon is classified as an "unusual spectrum 
showing a flat continuum with unusual features" (Chen et 
al, 2000). Aspin et al. (1988) and Yamashita et al. (1989) 
claimed to detect an extended emission eastward of R Mon, 
but it remained undetected by Close et al. (1997) with IRTF 
at 3.16 /im down to lcr ~ 0.05 mJy arcsec -2 (CVF filter 
with an exposure of 60 s). 

3 Spitzer observations 

We have obtained 4.2 hrs of Spitzer time to study the R 
Mon/HH 39 system. Our program (PID 20034) uses all in- 
struments: IRAC and MIPS images of HH 39; IRAC images 
of R Mon and its jets; IRS spectra of R Mon, its jets, and 
several of the HH 39 knots; and MIPS SED of R Mon. The 
observations were taken at different epochs, spanning from 
October 2005 to April 2006. 

3.1 The Herbig-Haro object HH 39 

Fig- HI (left) shows the contour profiles of the Herbig-Haro 
object HH 39 see in the R band (from Walsh and Malin, 




Fig. 1 (Left) The HH 39 Herbig-Haro knots in the R band (from 
Walsh and Malin, 1985). (Right) 3-color IRAC image of HH 39 (red = 
3.6 /im, green = 4.5 /im, blue = 5.8 /im). 



1985). Several knots are labeled, and in particular knot A has 
been identified as the working surface of the jet onto the sur- 
rounding molecular material. The right panel of fig. Q] shows 
a 3-color IRAC image of HH 39 (red = 3.6 /im, green = 
4.5 /im, blue = 5.8 /im). Several knots are clearly detected in 
the mid-infrared. Fig. [2] shows false-color images of the HH 
39 knots in each of the four IRAC bands. Knots D, G, C, and 
E are clearly detected in all IRAC bands, and knot H is tenta- 
tively detected. A diffuse emission is observed near the po- 
sition of knot A in IRAC1 only. The feature could be due to 
H2 lines or 3.3 /im PAH. Since this feature is not detected in 
IRAC2, this is a "PAH"-free band, it is very likely that PAHs 
contribute to the diffuse emission near knot A. H2 lines, on 
the other hand, most likely contribute to the emission in 
knots D, G, C, and E. Supporting evidence comes from a 
faint 0-0 S(l) rotational H2 line detected in the IRS spectra 
at 17.03 ^m of the bright knots (0.7 - 1.0 x 10~ 20 W cm~ 2 ). 
In contrast, IRS spectra of knots A+A' do not show evidence 
of H2 line emission in excess of the nearby continuum. Note 
that faint [Ne II] emission at 12.81 /im is measured in the 
bright knots and at the position of knots A+A' (of the or- 
der of 0.5 x 10~ 20 W cm -2 ). At 8 /im, the bright knots are 
barely visible on top of a diffuse emission that sweeps in 
through the knots. This extended emission is likely of due 
to PAHs (7.7 /im) and is also seen in the dust continuum 
24 /im MIPS image (Fig. [3] and [4]). Its origin is unclear but it 
could be the upper wall of the NGC 2261 reflection nebula 
cavity. 

Fig. [3] shows the MIPS images of the HH 39 region at 
24 and 70 /im together with an optical Deep Sky Survey 
image. Dust continuum emission is detected across the HH 
39 knots, in particular at 24 /im. Nevertheless faint HH 39 
knots are also detected at 24 /im (Fig. [4]), but there is no clear 
evidence of the knots in the 70 /im image. On the other hand, 
it should be reminded that the FWHM PRF at 70 fim is 18", 
i.e., about the distance from knot A to the group of knots 
(G, D, E, C, H). Therefore, since the MIPS 70 jim image 
does not, apparently, observe the same diffuse emission as 
at 24 /im or 8 /im, it is possible that the faint emission at the 
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Fig. 2 (a) IRAC false-color images of HH 39. Background surface flux 
densities are of the order 0.30,0.35,2.1, 10.5 MJy sr" 1 in IRAC 1, 2, 
3, and 4, respectively (zodi estimates: 0.12,0.37, 1.72, 10.4 MJy sr -1 ; 
ISM light estimates: 0.15,0.20,0.56,4.79 MJy sr" 1 ). In con- 
trast, the bright knots have peak surface flux densities of about 
0.95, 1.1,3.2, 11.8 MJy sr" 1 . 

position of the cluster of knots is in fact partially due to [O 
I] 63 jim line emission from the knots. 

Emission from the upper part of NGC 2261 is also de- 
tected in MIPS at 24 jim. The bright emission SW of HH 39 
is the IR source IRAS 06362+0853. There is no evidence 
for diffuse emission in 2MASS around IRAS 06362+0853, 
suggesting that it is of interstellar cirrus nature. Surface flux 
densities at 24 jim are around 42 MJy sr -1 for the "blue" 
color regions and go up to 43.1 MJy sr -1 across HH 39. 
Peak surface flux densities for IRAS 06362+0853 are about 
44.4 MJy sr -1 and 46.7 MJy sr -1 for the dust emission 
of NGC 2261. At 70 jim, "blue" regions are about 31 - 
33 MJy sr -1 , whereas the emission across HH 39 is 48 — 
50 MJy sr- 1 , and 80 - 85 MJy sr" 1 near IRAS 06362+0853. 
The SPOT background estimates are 37.6 MJy sr -1 (zodi) 
and 2.96 MJy sr -1 (ISM) at 24 jim and 10.8 MJy sr -1 (zodi), 
21.9 MJy sr" 1 (ISM), and 0.2 MJy sr" 1 (cosmic) at 70 jim. 
The surface flux densities in the "blue" regions are of similar 
values as the estimated total. 

3.2RMon 

Fig. \5\ shows the IRAC images of the Herbig Be star R Mon 
and its immediate surroundings. Note that IRAC2 could not 
be used, even in sub-array mode, due to R Mon's bright- 
ness at 4.5 jim. The main goal of these observations was 
to detect any faint emission from the jets (NS direction) or 
from the circumstellar disk (EW direction). In particular, we 
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Fig. 3 MIPS 24 and 70 jim images (left and middle) of the HH 39 re- 
gion together with a DSS image (right) of R Mon, the reflection nebula 
NGC 2261 and HH 39, on the same scale. 




Fig. 4 IRAC 3.6 jim and 8 jim, and MIPS 24 jim and 70 jim false- 
color images together with the positions of HH 39 knots detected in 
IRAC1. The IRAC 8 jim and MIPS 24 jim detect PAH and dust emis- 
sion across HH 39, likely due to the upper wall of the NGC 2261 reflec- 
tion nebula cavity. In contrast, at 70 /im, the upper wall is not detected, 
and the HH 39 knots are not clearly detected (see text). 

aimed to determine whether a faint eastward extended emis- 
sion feature detected by Aspin et al. (1988) and Yamashita 
et al. (1989), but undetected by Close et al. (1997), could in- 
deed be detected with the highly sensitive IRAC detectors. 

The total (zodi and cirrus) background flux densities are 
0.35, 2.86, 15.5 MJy sr" 1 (based on SPOT). However, the 
brightness of R Mon is such that the PRF illuminates the full 
sub-array detector. Aperture photometry (with a radius of 10 
pixels = 12", requiring no aperture correction), the total flux 
densities for R Mon are 17.59 ±0.29 Jy, 30.41 ±0.40 Jy, 
and 31.37 ±0.36 Jy in IRAC 1, 3, and 4, respectively. The 
estimated background fluxes over a circle of 10 pixel ra- 
dius are negligible (0.004, 0.03, 0.17 Jy) compared to R 
Mon's brightness and the RMS uncertainties. Note that the 
above fluxes include any contribution from R Mon B, sepa- 
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Fig. 5 IRAC 3.6 {im, 5.8 /im, and 8 (Am images (log scale) of the 
Herbig Be star R Mon. IRAC2 could not be obtained, due to R Mon's 
brightness at 4.5 /im, even with a sub-array mode. 



rated from R Mon by 0.7" (Close et al., 1997). However, in 
comparison with R Mon, the companion is expected to con- 
tribute negligibly (0.0013,0.0085,0.035 Jy in JHK'; Close 
et al, 1997). The next step will be to subtract the PRF of R 
Mon to possibly detect emission from the disk and the jets. 
We will need to create IRAC sub-array PRFs from sub-array 
data of stars observed, e.g., in the FEPS program. 

Finally, Fig. [6] shows the IRS SH spectrum of R Mon 
(top) and its raw SED from 1 to 100 /im (bottom). For the 
latter, we used IRAC fluxes, the IRS SH spectrum, and the 
MIPS SED and complemented them with values from the 
literature (Close et al, 1997; and MSX). The IRS SH spec- 
trum shows a clear PAH feature at 11.3 /im and other faint 
features might be present as well. A detailed analysis is on- 
going to remove instrumental effects (e.g., defringing). 
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